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ABSTRACT 

We investigate the enrichment of the pre-solar cloud core with short lived radionuclides (SLRs), 
especially ^^Al. The homogeneity and the surprisingly small spread in the ratio ^^Al/^^Al observed 
in the overwhelming majority of calcium-aluminium-rich inclusions (CAIs) in a vast variety of primi- 
tive chondritic meteorites places strong constraints on the formation of the the solar system. Freshly 
synthesized radioactive ^^Alhas to be included and well mixed within 20kyr. After discussing var- 
ious scenarios including X-winds, AGB stars and Wolf-Rayet stars, we come to the conclusion that 
triggering the collapse of a cold cloud core by a nearby supernova is the most promising scenario. 
We then narrow down the vast parameter space by considering the pre-explosion survivability of such 
a clump as well as the cross-section necessary for sufficient enrichment. We employ numerical sim- 
ulations to address the mixing of the radioactively enriched SN gas with the pre-existing gas and 
the forced collapse within 20kyr. We show that a cold clump of 10 Mq at a distance of 5pc can be 
sufficiently enriched in Aland triggered into collapse fast enough - within 18kyr after encountering 
the supernova shock - for a range of different metallicities and progenitor masses, even if the enriched 
material is assumed to be distributed homogeneously in the entire supernova bubble. In summary, we 
envision an environment for the birth place of the Solar System 4.567 Gyr ago similar to the situation 
of the pillars in M16 nowadays, where molecular cloud cores adjacent to an HII region will be hit by 
a supernova explosion in the future. We show that the triggered collapse and formation of the Solar 
System as well as the required enrichment with radioactive ^^Al are possible in this scenario. 
Subject headings: Stars: formation, Stars: protostars, (Stars:) supernovae: general. Hydrodynamics, 
ISM: kinematics and dynamics, ISM: abundances, Protoplanetary disks. Meteorites 



1. INTRODUCTION 

The time-scale for the formation events of our Solar 
System can be derived from the decay products of ra- 
dioactive elements found in meteorites. Short lived ra- 
dionucfides (SLRs) such as ^SAl, ^iCa, '^^Mn and *^°Fe 
can be employed as high-precision and high-resolution 
chronometers due to their short half-lives. These SLRs 
are found in a wide variety of Solar System materials, 
including calcium- alumini um-rich inclusions (CAIs) in 
prim itive chondrites (e.g. iLee et aDll976t lAmelin et al.l 
1200^. The decay of ^^Al into ^^Mg in particular, with 
a half-life of r « 0.7 Myr, provides by far the highest 
resolution measurements that m ark the initial forma - 
tion of the proto-planetary disk (jJacobsen et al.l 120081 ). 
As ^''Al decays, the ratio of ^'^Al/^'^Alat the time of 
condensation can be directly derived from the present 
day ^^Mg/^*Mg and ^''Al/^^Mg measurements. Further- 
more, the spread of ^^Al/^''A1 among different CAIs 
gives the precise time-span in which these CAIs formed. 
While there is considerable spread among the CAIs to- 
wards lower ratios of ^^Al/^^Al due to remelting and 
thermal metamor phism, the upper value, the so called 
'canonical ratio'. (iMacPherson et aLlll995D is now well 
estab lished (e.g. iJacobsen et al.l 120081 : IVilleneuve et al.l 
[200I. 

The general picture we adopt here is that a certain 
amount of ^^Al is injected in the nascent solar neb- 
ula and then gets incorporated into the earliest formed 



CAIs as soon as the temperature drops below the con- 
densation temperature of CAI minerals. These CAIs are 
most frequently found in CV-chondrites^. Lower ratios of 
^^Al/^^Al can then be explained by subsequent episodes 
of CAI remelting or thermal disturbance, thereby ex- 
plaining the heterogeneity below the canonical value. 
Therefore, the CAIs found in chondrites represent the 
first known solid objects that crystalized within our So- 
lar System and can be used as an anchor point to deter- 
mine the formation time-scale of our Solar System. Var- 
ious measurements of different CAIs by several research 
groups have not only confirmed the canonical ratio of 
(5.23 ± 0.13) X 10~^, but also e stablished a very sma, ll 
spread (for a recent review see iVilleneuve et ani2009t ). 
This spread corresponds to an ag e range of less than 
~ 20kyr (e.g. IJacobsen et al.|[2008f ). 

In addition, the Mg-isotope composition appears to be 
fairly uniform among bulk chondrites as well as Mars, 
Moon and Earth. This is a strong indication that initially 
Al- and Mg-isotopes are distributed quite homogeneously 
(±10%) in the p roto-planet ary di sk IjThrane et alll2006l : 
IVilleneuve et al . 2009; Bossi 2011[ ). The extremely smafi 
time-span together with the highly homogeneous mixing 
of isotopes poses a severe challenge for theoretical models 
on the formation of our Solar System 

^ CV-chondrites are a class of CA I-rich carbonaceous chondrites 
named after the Vigarano meteorite lIDauphas fc Chaussidonl2011l ) 
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Various theoretical scenarios for the formation of the 
Solar System have been discussed. Shortly after the dis- 
covery of SLRs, it was proposed that they were injected 
by a nearby massive st ar. This can happe n either via 
a supernova explosion (Cameron & Truran 1977') or by 
the s trong winds of a Wolf-Rayet star ( Arno uld et al.l 
Il997f ). Another possibility would be the in situ enrich- 
ment inside the disk b y the non-thermal a. c tivity of the 
Sun (jShu et all [20010 . iWasserburg et al.l ()1995[ ) have 
proposed an asymptotic giant branch (AGB) star as an 
possible source of SLR s in the Solar System (see also 
iTrigo-Rodriguez et aLllTooaV 

Up to now, none of these models have been able to 
explain the extremely small range of CAT condensation 
ages (the spread of w 20kyr). In this paper, we discuss 
several possible solutions to this problem. First, we show 
that most of these scenarios are not compatible with an 
almost coeval formation of CAIs. We then go on to per- 
form numerical simulations on the promising scenario of 
a type Ila supernova as the trigger of our Solar System. 

2. PHYSICAL PROCESSES AND FEASIBILITY OF 
SOLUTIONS 

In this section we discuss the physical processes in dif- 
ferent scenarios and their viability. For a detailed discus- 
sion of these scenarios, including an assessment of their 
likelihood see the recent review bv I Adams! dMol) . 

As mentioned before, the canonical value of ^^A1/^''A1 
is established to be (5.23±0.13) x 10~^. Assumi ng an to- 
tal mass for 27 Al of M27A1 = 7.1x10-5 Mq (e.g. lLoddersI 
[200l we can convert the observed ratio into a total mass 
of M^^I'Xr'''''' ~ 3.71 X 10-9 Mo of radioactive ^^Al, which 
has to be present in the early Solar System. 



2.1. In situ Enrichment 

One solution t o the inco r porat ion of ^^Al into CAIs 
was proposed bv IShu et al.l (|2001[ ). Young stars can un- 
dergo the so-called X-wind phase. In this phase the com- 
bination of magnetic fields, inflowing material and stellar 
outflows can form a w ind bet ween the proto-star and the 
proto-planetary disk (jShu et al. 1997). 

SRLs are produced in this scenario by solar ener- 
getic particle irradiat ion of dust and gas near the proto- 
Sun (|Shu et al.ll2001[ ) and incorporated into CAIs. This 
mechanism can explain the ratio of ^°Be/^Be = 8.8 ± 
0.6 X IQ-'' inferred fr om measurements of ^ °B, the de- 
cay product of ^°Be ()McKeegan et al.ll2000l). Another 
way t o explain this ratio was proposed by Desch et"all 
()2004|) . They calculate the contribution of ^°Be trapped 
by galactic cosmic rays to the collapsing molecular cloud 
core and conclude that a large amount, if not all, ^^Be 
can be attributed to cosmic rays. 

For the purpose of this study it is sufHcient to note that 
this mechanism can not explain the observed amount of 
2^A1 due to the following reasons. First o f all, Aland 
^°Be do not correlate well in CAIs (e.g. iMarhas et al.l 
|2002|) . In addition, SLRs produced by this mechanism 
are intrinsically expected to be highly heterogeneous. 
Variations in their relative abundances would reflect the 
local energetic particle environment and their episodic 
nature of production. This is in contradiction to the ob- 
served homogeneity (±10%) of ^^Al. Furthermore, the 
amount of ^^Al which can be produced by this mecha- 
nism is at least an order of magnitude too low to explain 



the canonical value (|Duprat fc Tatischefj|2007ll2008[ ). 

2.2. External Enrichment by AGB stars 

Another sou rce of Al ar e AGB stars. Every star of 
~ 0.8 - 8Mq (H erwi gl [2005h under goes this phase when 
helium fusion in the core is complete and helium shell 
burning begins. In this phase, the star loses most of 
its envelope in stellar winds, which could in te rm en- 
rich the proto-solar cloud(e.g. IWasserburg et al.l [19951. 
This stage is at the end of a stellar life. However, AGB 
stars are not present in star formin g regions and henc e 
an encounter is very unlikely (Kas tner fc Myer"slll994l l. 
In addition, the total amount of ^^Al requires a massive 
AGB-star and the short enrichment-timescale requires a 
brief AGB-phase. 

2.3. External Enrichment by Wolf-Rayet stars 

An interest ing possibility was proposed by 

iTatischeff elaH (1201011 wh o followed up on an ear- 
lier idea by lArnould et al.l ()1997l ). Here, a run-away 
Wolf-Rayet star that gets ejected from its parental 
cluster of massive stars sweeps up a smaller shell while 
traveling into the ISM. This shell gets enriched by 
material from the strong stellar winds of the Wolf-Rayet 
star. As soon as this star explodes in a supernova it 
sweeps up an even denser shell which expands rapidly. 
After ~ lOkyr the shell starts to cool and a cold core 
can start to form due to dynamical and gravitational 
instabilities which subsequently undergoes collapse. 
During this collapse phase the the surface density of 
the core increases with time (jWhitworth et al.l 11994 
iHeitsch et al.ll2008l ). Since this theory assumes that the 
2^Al/27Alratio increases as the surface density of the 
clump increases, the core becomes enriched sufficiently 
to explain the observed values in CAIs after 10^ — 10^ yr. 
However, there is no physical reason why this ratio 
should increase as the core density rises. The initial 
ratio within the shell should be conserved and only 
modified b y radioactive decay (t he decay is taken into 
account bv iTatischeff et al"1 l2010f). 

A similar scenario was proposed by iGaidos et al.l 
11009). Here, an entire molecular cloud (M > 10^ M©) is 
enriched in ^6 Al by a nearby star forming region. After 
^4 — 5 Myr this cloud begins to form stars, including 
the Sun. During this stage, the Solar System is enriched 
with heavier SLRs. This scenario can explain the differ- 
ent abundances of SLRs (see ijl]). However, it requires a 
precise chain of events. In addition, this theory assumes 
continuous enrichment which is at odds with the small 
age spread of CAIs (see below). 

2.4. Constraints on continuous enrichment 

All three scenarios described above are based on the 
assumption of a continuous enrichment. In that case, 
the characteristic time scale for the enrichment with ra- 
dioactive isotopes is much longer than the brief duration 
(< 20kyr) of the CAI formation episode which is inferred 
from the meteoritic data. 

In external enrichment models, sufficient radioactive 
isotopes can only be intercepted by a tenuous and ex- 
tended solar-system progenitor cloud prior to its collapse. 
But the free fall time scale for such a cloud is generally 
^ 100 kyr or longer. If the central star and disk ac- 
crete material from this parent envelope on such a time 
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scale, dust and grains would be processed with different 
rates at different locations and in different stages. Dur- 
ing this entire phase, the assembly and formation of the 
CAIs would be possible, leading to inhomogeneities and 
an age-spread of ~ fOOkyr. 

This is inconsistent with the observationally inferred 
homogeneity of Mg- and Al-isotopes during the forma- 
tion epoch of CAIs. Here, we make a clear distinction 
and differentiate the initial formation epoch of CAIs from 
later epochs when some CAIs subsequently experienced 
a complex thermal history within the solar nebula and 
metamorphism on their host m eteorite parent asteroids 
extending over 0. 1 - SMyr (e.g iMacPhe rson et al1ll995l: 
Hsu et al. | 120001 : lYoung et all [2005: M acPherson et al.l 



2007f) . One possible mechanism to homogenize the iso- 



tope's distribution and reset the clock for grain con- 
densation and CAI formation is to heat the infalling 
gas above the conden sation teniperat u re for CAIs o f 
~ 1500 - 1800K (e.g. iGrossmanI l200l iLoddersI [200l . 
During the infall, accretion and shocks can dissipate suf- 
ficient energy and heat up the disk to this temperature 
close to the central star (within ~ 0.1 AU). In order for 
the infalling gas to arrive at this location and form a 
compact disk, it must carry very little angular momen- 
tum. However, turbulent angular momentum transport, 
induced either by gravitational or magneto-rotational in- 
stabilities, would lead to mass diffusion, ^ yhen the disk 
spreads beyond an AU, its mid-plane cools (jRuden fc LinI 
Il986t iGaraud fc LinI l2007t ) well below the condensation 
temperature of the CAIs. If the infalling gas has a mod- 
est amount of angular momentum, it would form more 
extended disks with shoc k temperatures bel ow lO'^K at 
a distance outside 1 AU (jWalch et al.l 120091 ). In either 
limit, condensation of the elements and the formation of 
CAI's would continue until infall is terminated. 

A possible solution to reconcile with CAI's brief forma- 
tion epoch is to assume that all infalling gas was heated 
above their condensation temperature of ~ 1800K and 
then cooled within 20 kyr. Such a scenario would require 
a very compact disk and therefore a violent, external 
cause, i.e. a supernova, to trigger and drive the collapse 
of the entire progenitor cloud which formed the Solar Sys- 
tem. As shown in Appendix A, the subsequent growth 
of CAIs would be possible on a short enough time-scale 
in such a disk. 

2.5. External Enrichment by a Type 11 Supernova 

A frequently discussed possibility for SLRs is the ex- 
ternal enri c hmen t by a type Ha supernova (SN Ha). 
iHuss et al] ()2009t ) estimate that a supernova with a pro- 
genitor of 20 — 60 Mq can produce the required SLRs. 
From the simple assumption that the cross-section of 
the enriched region has to be high enough to intercept 
enough ^^Al we can derive a size-distance relation from 
geometrical considerations: 
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MA126 



Ma 126 



(1) 



Where R is the radius of the intercepting cloud or disk 
and D is the distance from the supernova. M^jgi and 

Af^i^e^ are the amount of ^^Al required in the Solar Sys- 
tem and produced in the supernova, respectively. This 
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Fig. 1. — The parameter space. Plotted is the distance D from 
the exploding star versus the radius R of the pre-solar clump. The 
diagonal lines give the required ratio for the enrichment. Black is 
for a 40 Mq star, blue for a 20 Mq star, solid is solar metallicity 
with Z=0.02, dashed is Z=0.004. The vertical lines are the radius 
of the Hll-region for an 05V star after SMyr and a BOV star after 
lOMyr, respectively - dotted is for a cold medium of 5 X 10* cm~^, 
dash-dotted for 5 X 10'^ cm~^. The diamonds represent the param- 
eters chosen for the simulations in ijS] 



M 



Mass [M^ 



A i26 



'M 



0J 



Z=0.02 



Z=0.004 



20 
40 



1.5 X 10-5 8.65 X 10-6 
6.64 X 10-5 2.23 x 10-^ 



TABLE 1 

Typical ^^Al yields in Mq for different supernova 

PROGENITOR MASSES AND META LLICITIES TAKEN FROM 
INOMOTO ET AL.I ll2006i) . 



equation is of course simplified, e.g. neglecting radioac- 
tive decay as the time-scales are very short, as well as 
assuming a perfect mixing efficiency, i.e. all ^^Al inter- 
cepted is assumed to be mixed into the Solar System. 



2.5.1. Enrichment in the Disk Stage 

From the values in Table |T] one can calculate that at 
solar metallicity a 40 Mq progenitor supernova could de- 
posit enough ^^Al in a 1000 AU face-on disk at D — 
0.3 pc or a 20 M© progenitor enough ^^Al at D = 0.15pc. 
That is very close to the hot, ionizing star. Although a 
proto- stcUar disk could survive in this violent environ- 
ment (lAdamj|2010i). there are additional timing issues. 
For example iGounelle fc MeibomI (|2008l ) estimate the 
probability that the proto-solar system is young enough 
(< IMyr), able to survive photo-evaporation and close 
enough to get sufficiently enriched {D < 0.3 pc) to be 
less than 10"'^ in the most favorable case. Based on the 
requir ed young age of the system, IWilliams fc GaidosI 
(|2007( ) argue for a very massive supernova progenitor, 
which in turn leads to a short life-time, i.e. an early su- 
pernova. As very massive stars are only found in massive 
cl usters, they con c lude that th e likeli hood is < 1%. 

lOuellette et"all (|2005L I2009L I2010D suggested an in- 
triguing possibility where the disk is further away from 
the supernova ('^ 2 pc) and the enrichment is achieved 
via the injection of already formed grains into the pre- 
solar disk ('aerogel' model). To achieve the high ob- 
served yield, the supernova gas has to be highly clumpy 
and the clumps have to be highly enriched in SLRs - 
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to this possibility the authors estimate a probabihty of 
^ lO^'^ — 10^^. There are, however, additio nal uncertain- 
ties. Recent models fe.g. Ilda fc Liij[200l) suggest that 
the minimum mass solar nebula might be more massive 
than their assumed 0.01 Mq, which would require an even 
more efficient enrichment. Future investigation on the 
dumpiness of SN-shock fronts are desirable (see A 
remaining issue in this picture is the question why this 
scenario does not result in multiple enrichments of the 
early Solar System. After the first enrichment by a su- 
pernova, which leads to the canonical value of ^^Al/^^Al, 
it is likely that several more lighter and therefore more 
long-lived stars in the cluster undergo subsequent super- 
nova. This would result in several enrichment episodes 
during the entire lifetime of the disk (> lOMyr). 

2.5.2. Enrichment in the Cloud Core Stage 

The enrichment of the pre-solar cloud core by 
a supernova has b een s uggested for a long time 
(ICameron fc TruranI 119771: ICameron et al.l Il995l[ . Re- 
centlv.'B oss et al.l ()2008ll201 0) and B oss fc Keiseij (|2010[ ) 
investigated the implosion and enrichment of a 1 Mq 
cloud core in detail. They assume the cloud core to 
be more than 10 pc away from the progenitor. There- 
fore, the supernova front has slowed down considerably 
by the time it reached the pre-solar cloud. In addi- 
tion, it has a finite thickness. In their models, the shock 
speed ranges fromws = l-100kms~^ ()Boss et al.ll2010[ ). 
the shock thickness from = 3 x 10"'* — 3 x 10~^ pc 
and the shoc k density from ps = 3.6 x 10^^^ — 2.4 x 
lO-i'^gcm-3 (|Boss fc Keiseij|2010[ ). They find that with 
certain initial conditions the enrichment is possible in 
their model. For example a thin (Js = 3 x 10"'*), dense( 
ps — 6 X 10~**gcm~^) shock at a rather low velocity 
(us — 20kms~*) and therefore far away from the source 
can enrich the clump sufficiently. This can be easily un- 
derstood as they assume that the enriched supernova 
material is completely within the shock. Therefore, as 
the geometrical dilution rises, the shock has to be thin- 
ner to contain a sufficient amount of supernova mate- 
rial. In addition, the time-scale of the collapse is longer 
{t > 40kyr) than the value inferred from CAIs in this 
scenario. Again, this can be attributed to the slower su- 
pernova shock further away from the progenitor. Never- 
theless, this is an encouraging development and we think 
this scenario is the most promising. 

Bearing the geometrical dilution in mind (cf Fig. [T]) 
we revisit this scenario but place the pre-solar clump 
closer to the supernova. Therefore, the shock velocities 
are much higher (see §3.2.1|) . We also take the more 
conservative assumption that the material is spread out 
in the hot phase or at least in a ring at the border within 
the supernova explosion bubble (see ii3.2.2l and I3.2.3p . 
Furthermore, we do not assume the shock already formed 
a density enhancement. Instead, we setup the supernova 
blast wave in the velocities under the assumption of a 
thin shock. As the evolution is expected to be much more 
violent, we investigate the collapse of a heavier (10 M©) 
cloud core (see i i3.2.ip . To assess the feasibility of this 
initial condition, we address the aspect of survival of the 
pre-solar cloud ( §3.ip . 

A slightly different approach would be the enrichment 
in the cloud stage without triggering the collapse. Then, 



the cold cloud would collapse in an 'isolated' manner at 
a later stage. However, there are clear short-comings 
of this scenario. To allow for the tight spread of only 
~ 20kyr, the collapse would have to lead to a stage, 
where the entire disk is heated above ~ 1800 K and then 
homogeneously cools to lower temperatures within 20 kyr 
(see also ij2.4p . To achieve the required heating from 
the gravitational energy the accretion disk would have 
to be < 0.5 AU. At this scale, the disk could then cool 
coherently, the CAI could condense and later migrate 
outwards to explain the homogeneous ^^Al/^^Al in CAIs 
in different meteorites. 

However, isolated collapse leads to larger, cooler disks. 
Thus, a triggering event will be required to force the col- 
lapse to a small scale. It is straightforward to assume 
that this triggering event was the same event as the en- 
richment, i.e. the supernova shock. Thus, we investigate 
the triggered case in more detail. 

3. MODEL AND SIMULATIONS 

In this section, we address the parameter space by ac- 
cessing the enrichment, survivability and existence of the 
pre-solar cloud core in or adjacent to the HII region of the 
massive star. We investigate the mixing and the collapse 
probabilities by employing numerical simulations. 

In Fig. [T]we plot the distance of the clump to the pro- 
genitor star versus the size of the clump. From the ge- 
ometrical enrichment cross-section discussed before, we 
can derive a first constraint. Taking the yields in Table 
[1] into account we plot Eq. [TJ These diagonal lines de- 
note the maximum enrichment efficiency. Even if all the 
material intersected by the initial cloud ends up in the 
Solar System (100% efficiency), the initial condition has 
to be above these lines to explain the ratio inferred from 
CAIs. 

3.1. Survival 

The most obvious proof for the existence of cold cloud 
cores in the proximity of massive sta r comes from obser- 
vations. In the pillars in M16 ( e.g. iHester et al.lll99^ 
iMcCaughrean &: Anderseiil (|2002l ) find several cold cores 
in almost every region of the pillars. These pillars are 
about 2 pc away from the ionizing so urces, the age of th e 
region is estimated to be 1 - 2Myr (jHester et al.lll996l l. 
This agrees well with our assumption of structures being 
present at 5pc after 3Myr (see below). 

From a theoretical point of view, recent simulations 
on the formation and evolution of Hll-regions show 
the u biquity of these struc tures (e.g. [M cUcm a et al.l 
2006'; 'Krumholz et aL| l2007t iGritschneder et al.i l2009bl: 
Arthur et al. 2011). These simulations focus mostly on 
th e ionization of turbulent molecular clouds. As shown 
in IGritschneder et al.l ()2010[ ). the average time evolution 
of the front position is sti ll wel l appr oximated by an an- 
alytic solution as given in ISh^(ll99l : 

4 

i?(i)=i?s(l + J^(i-<o))\ (2) 

where Rg is the initial Stroemgren radius^ and Cg^hot is 
the sound speed of the hot, ionized gas. Using this equa- 

^ The Stroemgren radius is the radius which can be immediately 
ionized by an ionizing source, without taking any hydrodynamic 
evolution into account l IStromgrerJIlQSSf) . 
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tion, we can estimate the size of the Hll-region. To inves- 
tigate the most extreme cases, we assume a short-hved 
40 M0-star with a very short hfe-time of 3 Myr as well as 
a 20 M f7)-star with a life-time of 10 Myr fe.g. lHurlev et al] 
12000^ in different density regions. The resulting radii are 
depicted by the vertical lines in Fig. [1] Surviving clumps 
should be to the right of theses lines. Note that we as- 
sume a relatively high density in the neutral medium 
before the ignition of the star in order to re flect the 
birth environment of massive stars. 'Mvers ( 2009) suggest 
that the birth of stellar groups is associated with hubs 
of column densities greater than 10^^ cm~^ , for example 
in Taurus ([Goldsmith et al.ll2008l ) , which corresponds to 
the assumed densities here. In addition, the high densi- 
ties have to be present only in the direct surrounding of 
the massive stars, as the further evolution given in Eq. [5] 
only depends on the initial Stroemgren radius (and the 
sound speed in the hot gas). Altogether, these values are 
chosen to reflect our current understanding of massive 
star formation. In addition. Fig. 1 gives the averaged 
front position, the pillars in the observations as well as 
the simulations can be a few pc further inside, pertruding 
into the ionized bubble. 

3.2. Collapse Timescale and Mixing in the Gas Phase 

The remaining issues are whether the gas from the 
supernova can be efficiently mixed with the gas from 
the pre-existing core, whether collapse occurs and if so, 
whether it is fast enough. To test these questions, we 
employ numerical simulations. 

3.2.1. Initial Conditions 

Because the main focus lies on mixing as well as cool- 
ing, we employ the numerical code COSMOS. The non- 
relativistic hydrodyn amics of COSMOS are discussed in 
(|Anninos et al.ll2003D . For the models described below, 
the code is run with the default hydrodynamics solver, 
which is a total variation diminishing (TVD) approx- 
imate Riemann solver a nd a t hird-order Runge-Kutta 
time- marching scheme C Shu fc O sher 1988). Because of 
the importance of cooling, the code is run in internal en- 
ergy mode. Shocks are handled with a zone-centered 
scalar artificial viscosity. Optically thin heating and 
cooling are included, using equilibrium cooling curves 
IPalgarno & McCray 1972; Bochringer & Henslcr 1989). 
Cooling by metals assumes solar metallicity. Mixing is 
followed by the use of a passive tracer, whose value is 
initialized to be proportional to the gas density behind 
the blast-wave. The problems are run in two dimensions. 
Cylindrical symmetry is imposed by taking advantage of 
the relativistic capability of COSMOS, and running the 
problems with a cylindrical metric. 

The problem region runs from -5 to 10 cloud radii in 
the direction of the shock, and from to 4 cloud radii in 
the transverse direction. The large radial size is adopted 
so as to prevent any boundary affects from altering the 
cloud evolution. Due to the motion of the cloud after 
shock passage, we use a uniform grid spacing in the di- 
rection of the incoming shock, with 4096 zones along the 
axis, or 274 zones across the cloud radius. In order to 
enhance resolution of the collapsed cloud, geometrically 
increasing zone spacing is used in the radial direction. 
A total of 512 zones are included radially, with a ratio 



of 1.005 in their spacing, giving a minimum zone size of 
0.0017 times the cloud radius. 

We set up an molecular cloud core in isolation, which is 
going to be hit by the supernova shock- wave. We approx- 
imate t he core by an marginally stable Bonnor-Ebert- 
Sphere ()Bonnoilll956l ). the solution of the Lane-Emden 
equation for hydrostatic equilibrium 



1 d dp 

— — — — = -AnGp, 
ar p ar 



(3) 



where gravity and hydrodynamical pressure balance each 
other. We assume an initial temperature of Tcom = 20K 
and a central number density of Uc = 5.7 x 10^ cm~^, 
which gives a cloud core of mass M = 10 M0 with a ra- 
dius Rq — 0.21 pc. The sphere is embedded in a warm gas 
in pressure equilibrium with the sphere at a temperature 
Twarm = 792 K. This is realistic, as the core is embedded 
in the interstellar medium close to an Hll-region. As this 
is an idealized case and the main focus of this study is the 
enrichment with SRL s, we don't take a ny rotational or 
turbulent motion fe.g. lWalch et al.l2010[ ) inside the cloud 
core into account. While this initial rotation might be 
important in the later phases of the collapse - especially 
in determining the final disk size - its influence will be 
very small in the initial phase here, as the rotational ve- 
locities are orders of magnitude smaller than the shock 
speeds. In addition, the total angular momentum carried 
by a minimum mass nebula (0.02 Mq) out to 30 AU cor- 
responds to that of 1 Mq of material at 0.6 AU. Unless 
there is a substantial angular momentum loss since the 
initial collapse to the setup of a minimum mass nebula, 
the initial cloud must carry very little angular momen- 
tum. In such a low-angular momentum cloud, rotation is 
unlikely to inhibit the collapse until the size of the cloud 
is well below the resolution of our simulations and the 
dynamical time scale is much shorter than the 20 kyr of 
interest. 

The supernova shock-wave at the distances consid- 
ered here (5 — 10 pc) has not yet reached the snow-plow 
phase. So the shock is still in the energy dominated phase 
and can be a pproxim ated by a Sedo y- Taylor blast wave 
( Tavlo j[l950t ISedov|[l959l) . Followine lmoslev fc Wea^ 
( 1995( 1 we assume an energy release of E = 10^"'^ erg 
for core-collapse supernovae with progenitor stars in the 
range of 20 — 40 Mq . Considering the parameter space 
(i jS.ip . we perform two simulations, placing the cloud 
core at a distance Z? = 5 pc (case I) and D = 10 pc (case 
II), respectively, from the supernova. We evaluate the Se- 
dov solution for the given distances. In case I, the results 
are that the blast wave reaches the core after ~ 8.5 kyr 
with a post-shock speed of Vs = 276kms~^. The sound- 
speed in the post-shock gas is Cg = 129 km s^^. In case 
II the cold cloud core is reached after ~ 48.2 kyr with 
Vs = 97.5 km s~^ and Cg = 45.4 km s~^. 

3.2.2. Case I 

In the first case, we place the cold pre-solar core at 
a distance I? = 5 pc from the supernova. Fig. [2] shows 
the time evolution of the density in this simulation. The 
shock wave is propagating from the bottom to the top. 
As it can be clearly seen, the shock wave encompasses 
the cold core rapidly. Various hydrodynamical insta- 
bilities occur at the interface. After t = 4.15 kyr the 
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Fig. 2. — The time evolution of case I. Color coded is the density at t - 
units of the radius of the initial cold core (i?o = 0.21 pc). 

front has already passed the core. The part of the front 
which hit the clump continues to interact with the pre- 
existing core, leading to two fronts wrapping around it. 
As these fronts collide, the material gets mixed with the 
pre-existing cold gas and is driven into collapse. After 
t = 8.33 kyr the central region is already at a very high 
density (pmax = 3 x lO^^^gcm^^), the computation gets 
expensive and we terminate it. 

Fig. [3] shows a zoom-in into the density at the final 
snapshot. In Fig. 2] we show the temperature distri- 
bution in this region. As one can already see, the cen- 
tral region is cold. We define the central region as the 
cold, dense material between 0.25 and 0.75 on the y- 
axis of Fig. |4l corresponding to an elongated filament of 
w 0.1 pc diameter. The mass in this region below a tem- 
perature of 20 K is Mcorc — O.ISMq, the average density 
is p = 1 X 10""'^'' g cm~'^. This gives a Jeans-mass^ of 
~ 0.08 Mq. The core is twice as heavy and therefore will 
definitely undergo collapse. This collapse will proceed at 
least with the free-fall time-scale, which is w 20 kyr at 
this average density. More likely, it will proceed faster as 
the central density is two orders of magnitude higher and 
since the material is already moving inwards. In sum- 
mary, it is probably safe to say that the further collapse 
to a stage below 1800 K, where the CAIs form, takes less 
than 10 kyr. 

The mass is still lower than the solar value, but given 
the fact that cooling and accretion will still happen this 
value can increase. Although, given from the mass avail- 
able in the surrounding we do not estimate it to reach 
more than 0.2 M© in total. In addition, very strong in- 
flows or colliding remnants of the shock are still present, 

^ The Jeans-mass is calculated with a mean molecular weight of 
fi = 2.5, i.e. assuming that molecular hydrogen and helium are the 
dominant species. 



: Okyr, t = 4.16 kyr and t = 8.33 kyr. The length scale is given in 

flowing towards the centre with velocities of ~ lOkms^^. 
These will lead to an enhanced, fast collapse likely force 
the material in a small central region, which favors the 
simultaneous crystallization of CAIs. 

The other main focus here is the mixing of the enriched 
hot supernova gas with the cold pre-existing clump. To 
follow this evolution, we include a tracer in the post- 
shock gas inside the shock-front. The ratio of Al to the 
total gas mass within the supernova bubble depends on 
several factors. First of all, different progenitor masses 
give different yields. Secondly, the exact distribution of 
the enriched gas within the supernova bubble remains an 
open question. 

We investigate two cases here. The worst case, where 
the supernova enrichment is fully mixed with the entire 
material in the hot bubble, i.e. the 20 M© of total su- 
pernova yield are diluted with ^ 1000 M© within the 
bubble^. The second case is more optimistic, assuming 
that the enriched material stays close to the blast-wave. 
We assume it to be within the outer 10% of the radius 
of the bubble, leading to a dilution by only 350 M© of 
pre-supernova material. 

However, even in the worst case of complete mixing, 
the enrichment is sufficient for a 40 M© progenitor star at 
solar metallicity. This is shown in Fig. [S] Color coded is 
the ratio of Alto total mass, normalized to the required 
value of ~ 3.71 x 10~^. Since in the collapsing regions the 
ratio is of order unity or higher, we conclude that the gas 
phase is sufficiently enriched to allow for the formation 
of the observed CAIs later on. 

In Table [£] we list the outcome of different scenarios. 
We give Mfg^p the total mass of ^^Al inside the final 

To estimate the dilution we assume an average density of n = 
100 cm~^ within the supernova bubble from here on. This is a high 
value, but still realistic (for a discussion see Q. 
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Fig. 3. — Zoom in into the final stage of Case I at t = 8.33 kyr. 
Color coded is the density, the contours denote the gravitational 
potential. The length scale is given in units of the radius of the 
initial cold core (Ro = 0.21 pc). As can be clearly seen the densest 
regions will undergo gravitational coUajjse. 

collapsing core, as well as the dilution /o (i.e. the fraction 
of supernova material injected into the final core) and the 
injection efficiency /j as defined in Eq. [B] Finally, we list 
the resulting enrichment for a collapsing core scaled to 
the solar ratio 



Mcore M: 



canonical ' 
5A1 



(4) 



/q can be interpreted as the success of the model, i.e. 
value of f.O corresponds to the observed ^^Al/^^Al in 
the Solar System. 

Our results differ depending on the assumed distribu- 
tion of enriched material inside the supernova bubble and 
on the progenitor mass. As can be seen, in several cases 
the enrichment is sufficient, if the cloud core is placed 
at D = 5 pc from the supernova initially. Note that 
the dilution /o is systematically lower than the values 
of /o « (0.5 — 6) X 10~^ as given by iTakigawa et alj 
()2008| ). At the same time the mass in the cold core is 
only 0.13 Mq in our case, explaining that difference of one 
order in magnitude. In addition, we assume a high num- 
ber density inside the supernova bubble (n ~ 100 cm~^), 
which leads to a very strong dilution. Thus, these esti- 
mates are very conservative and can be viewed as a worst 
case scenario enrichment. 



Fig. 4. — The temperature distribution at the final stage at t = 
8.33 kyr. The blue are cold regions, red are hot regions. The length 
scale is given in units of the radius of the initial cold core (Rg = 
0.21 pc). The dense regions are very cold, indicating gravitational 
collapse. 

To further understand the reason for this very efficient 
mixing, it is useful to take a closer look at the earlier 
stage of the simulation. In Fig. [5] we show the tempera- 
ture distribution at the intermediate stage (t = 1.25 kyr). 
As can be clearly seen the material is already cooling be- 
hind the shock front. This cold, enriched clumps are 
mixing very efficiently with the pre-existing core. There- 
fore, a precise prescription of cooling processes is vital to 
investigate enrichment scenarios. 

3.2.3. Case II 

In the second case, the cold cloud core is at a distance 
13 = lOpc and is embedded in an ambient medium which 
is a factor of ten less dense with a ten times higher tem- 
perature than case I. Therefore, the ambient pressure is 
the same. Due to the longer distance the angular cross- 
section is smaller, the enriched material is already more 
dispersed. In addition, the shock has already lost more 
energy, i.e. the front speed is slower. 

The effect of this can be seen in Fig. [T] After t = 
14.6 kyr the cold core is still much less deformed by the 
supernova and not driven into collapse. Although the 
mixing seems to be very efficient, a closer look at Fig. [5] 
shows that the gravitational collapsing region, a remnant 
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Distance 


Mass 


Metallicity SNe mixing 




©] 


/o 




/i 




/© 


5 c 
pc 


40 Mq 


Z = Q 02 


shell 


1.6 X 10" 


-S 


2.3 X 10" 


-5 


5.3 X 10" 


-2 


3.34 








complete 


4.4 X 10" 


■10 


6.6 X 10" 


-6 


1.5 X 10" 


-2 


0.91 






Zj — U.UU'i 


sneii 


5.4 X 10" 


■10 


2.3 X 10" 


-5 


5.3 X 10" 


-2 


1.1 








complete 


1.6 X 10" 


■10 


7.0 X 10" 


6 


1.6 X 10" 


-2 


0.32 




20 M0 


^ — Q Q2 


shell 


3.6 X 10" 


10 


2.5 X 10" 


■5 


5.6 X 10" 


-2 


0.75 








complete 


1.0 X 10" 


10 


6.6 X 10" 


6 


1.5 X 10" 


-2 


0.21 








sneii 


1.6 X 10" 


10 


2.5 X 10" 


-5 


5.6 X 10" 


-2 


0.43 








complete 


5.7 X 10" 


11 


6.5 X 10" 


-6 


1.5 X 10" 


-2 


0.11 


10 pc 


40 Mq 


^ — Q Q2 


shell 


3.9 X 10" 


12 


5.8 X 10" 


-« 


5.2 X 10" 


-4 


0.065 








complete 


1.1 X 10" 


12 


1.6 X 10" 


-8 


1.5 X 10" 


-4 


0.017 






Z = 0.004 


shell 


1.3 X 10" 


12 


5.8 X 10" 


-8 


5.2 X 10" 


-4 


0.022 








complete 


3.8 X 10" 


■13 


1.8 X 10" 


■8 


1.6 X 10" 


-4 


0.0064 




20 M0 


Z = 0.02 


shell 


8.7 X 10" 


13 


5.8 X 10" 


■8 


5.2 X 10" 


-4 


0.014 








complete 


2.4 X 10" 


13 


1.6 X 10" 


-8 


1.5 X 10" 


-4 


0.0041 






Z = 0.004 


shell 


5.1 X 10" 


13 


5.9 X 10" 


-8 


5.4 X 10" 


-4 


0.0085 








complete 


1.4 X 10" 


13 


1.6 X 10" 


-8 


1.4 X 10" 


-4 


0.0023 



TABLE 2 

Outcome of different scenari os. We give the dista nce, mass, chosen metallicity and assumed SNe mixing. The 

CORRESPONDING YIELDS ARE TAKEN FROm INOMOTO ET AL.I (120061 '). We GIVE M|g™j, THE TOTAL MASS OF -^^Al INSIDE THE FINAL CORE, AS 
WELL AS THE DILUTION /q AND THE INJECTION EFFICIENCY /j . FINALLY WE LIST THE RESULTING ENRICHMENT FOR A COLLAPSING CORE 
SCALED UP TO THE SOLAR MASS, I.E. A VALUE OF 1.0 OR BIGGER STANDS FOR SUCCESSFUL ENRICHMENT TO THE SOLAR VALUE. 



of the original cloud core is not strongly enriched (the 
blue region at 0.0 — 0.1 cloud radii showing the highest 
density in Fig [7]) , whereas the enriched region further 
away (the light blue region between 0.7— 1.3 cloud radii) 
shows the effect of the two flanks of enriched supernova 
material colliding. 

In addition, the mass in the cold region is much 
smaller, only 0.01 M0, which is 30 times less than the 
Jeans mass at the average density of p = 7 x 10~^^ g cm~'^ 
in this region. Therefore, the cloud will not collapse in 
this case. This is due to supernova shock wave already 
being substantially weakened so far away from the front. 
Of course it would be possible to start with a more con- 
densed core, but then the enrichment, which is already 
difficult in this case (see below), becomes even more un- 
likely as the radius of the core would be smaller (cf Eq. 

From Table [2] it can be seen that sufficient enrichment 
is very unlikely in this scenario, especially since the en- 
riched material is now more diluted, as the volume of the 
supernova bubble is bigger, i.e. the material produced by 
the supernova is now mixed with ^ 10^ Mq in the com- 
plete mixing case and with ~' 2800 Mq in the shell case. 

Altogether the clump will neither collapse, nor be suf- 
ficiently enriched in this case. 

3.3. Comparison of the enrichment in simulations with 
theoretical (geometrical) predictions 

Here, we compare the predicted enrichment from the 
geometrical cross-section (Eq. [1] and Fig. [IJ with the 
finally achieved enrichment. 

From a geometrical point of view, the dilution /o is 
given as by the radius of the initial sphere Rq and the 
distance D to the source as: 

JO.gcora - 4£)27^- 

Thus, /o,goom is the amount of ■^^Al intercepted geomet- 
rically by the initial Bonnor Ebert sphere. /o,gcom ~ 
4.4 X 10"'' and /o.gcom = 1.1 x lO""' for cases I and II, 
respectively. The values are systematically higher than 



/o in Table [2] as Eq. [5] assumes perfect mixing of the su- 
pernova material with the background gas. However, it 
can be used to determine an injection efficiency, which we 
can not define directly, as we do not assume the material 
to be in the front initially. 

Instead, we define the injection efficiency /i as the ratio 
of the finally incorporated mass to expected intercepted 
mass from geometrical considerations 



f fl/fSNc f ■ ^ ' 

7 0,gcom-i^-'26^j jO.gcom 

/i would be unity if all intercepted material is included 
in the collapsing core. We give the values in Table [2] as 
well. For case I /i « 0.01 — 0.05, i.e. a few percent of the 
material crossing the cold sphere initially are incorpo- 
rated in the collapsing core finally. In case II, the values 
are much lower, fi « 0.0001 — 0.0006. First, this is due 
to a larger amount of un-enriched material in the bigger 
bubble. Second, the shock front is much weaker in case 
II, and therefore the front does not mix as efficient with 
the cloud core. 

4. DISCUSSION 

The outcome of these simulations and the results given 
here are a first order proof of principle. From the current 
scenario it seems possible that every supernova with a 
progenitor of 20—40 Mq can produce a system with solar- 
like ^^Al abundances by triggering a nearby (« 5pc) 
cloud core into collapse. However, we only tested a 
small range of the parameter space. Future simulations 
are desirable to further probe the likelihood. A first 
issue would be that we assume the (idealized) analyti- 
cal Sedov-Solution for the blast wave. A more realistic 
treatment for the entire evolution of the shock from the 
progenitor is required in the future. This is of course 
computationally very expensive. Nevertheless, it is im- 
portant to test e.g. the fragmentation and cooling of the 
shock front during the 8 kyr before it hits the cloud core. 
Another area to improve would be to include more de- 
tailed effects. For example, the core could be initially 
rotating. We assume here that the high shock velocities 
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Fig. 5. — The tracer field, which follows the supernova gas at t = 
8.33 kyr. Color coded is the ratio ^^Al/^^Al, the required value as 
inferred from CAI measurements is unity (green). The length scale 
is given in units of the radius of the initial cold core (Ro = 0.21 pc). 
The color coding here is for a 40 Mq progenitor at 5 pc distance and 
the supernova gas getting fully mixed within the entire supernova 
bubble. The contours show the density distribution. As it can be 
seen the collapsing pre-solar gas gets sufficiently enriched in this 
case. 

dominate over the rotation, but this should be further 
investigated. Furthermore, we neglect magnetic fields. 
They are of course present and should be involved in 
various processes, e.g. the details of the shock front and 
especially later on in the disk formation and the remove- 
ment of angular momentum. In the early stages, how- 
ever, the kinetic energy of the shock plays the dominant 
role. Still, it is worthwhile to investigate this regime of 
the parameter space in future simulations. 

A remaining issue is the production of heavy SLRs 
like ^^Mn and *^°Fe. iGounelle et"all (|2009t) suggest that 
the observed amount of ^"Fe can be inherited from the 
ISM. Still, the amount of galactic ^^Al is an order o f 
magnitude lower than in the Solar System (|Diehl|[2"006[ ). 
Thus, a different enrichment mechanism for ^^Al is still 
needed. In addition, in the supernova-enrichment sce- 
nario, ^°Fe and especially ^'^Mn are over-produced by 
a factor of 10-100 relative to their abundance in me- 
teorites when compared to the lighter SLRs ^^C a and 
^^AUe.g. lRauscher et al.ll2002t iGaidos et al.ll2009[ ). This 
is often used as an argument against the supernova- 
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Fig. 6. — The temperature distribution at a much earlier stage 
{t = 1.25 kyr). Blue is cold, red is hot. The length scale is given in 
units of the radius of the initial cold core {Rq = 0.21 pc). As can 
be seen, the enriched gas cools behind/within the shock- front and 
only later gets mixed with the pre-existing cold gas. The mixing 
happens between the two cold phases, rendering it more efficient. 

enrichment scenario. However, there are several possible 
solutions. First of all, the precise determination of abun- 
dances, especially of radioactive isotopes in supernovae 
remains challenging. Second, there are alternative mod- 
els, e.g. fallback-mixing within supernovae. Here, the 
innermost layers (commonly within the Si-burning layer) 
fall back onto the star. Therefore, there are much less 
SLR s he avier than S i within the supernova shock (e.g. 
iTakigawa et al.ll2008 ). Third, as it is an overproduction, 
not a lack of these SLRs, there could be different incor- 
poration efficiencies into CAIs and meteorites, which are 
not yet well understood. 

On a side note, the distribution of the supernova yield 
in the entire bubble versus the distribution of it in an 
outer rim gives different results (cf Table [2]) . As the iron 
is produced in layers further inside it is valid to speculate 
that the different ratios might be due to an homogeneous 
distribution of ^°Fe in the entire bubble, whereas ^^Al 
might be concentrated in an outer layer. This scenario is 
of course only valid if the initial distribution of elements 
in the supernova is conserved while the bubble expands 
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Fig. 7. — Zoom in into the final stage of Case II at t = 14.66 kyr. 
Color coded is the density, the length scale is given in units of the 
radius of the initial cold core (Rg = 0.21 pc). Not even the densest 
regions undergo gravitational collapse. 

to 5 pc. 

The different abundances of SLRs are often used to 
calibrate tlie delay between the supernova and the for- 
mation of CAIs. Mostly, these calculations give val- 
ues of about a Myr (Takigawa ct al. 2008: 0.7 Myr, 
iLoonev et al] 120061 : 1.8 Myr). This contradicts the fast 
formation time required to produce the small spread ob- 
served in the CAIs. In our opinion, the supernova yields 
of SLRs are at present not understood well enough to 
derive meaningful delay time scales between supernova 
and CAI formation. Especially, since there might be dif- 
ferent processes during the disk formation, which could 
also explain the different abundances of SLRs. 

In our simulations we find that a distance 13 w 5 pc is 
optimal to keep the balance of survivability for the pre- 
existing cold clump and the sufficien t enrichment. This is 
larger than previous estimates - e.g. ILoonev et al.l (|2006[ ) 
derive a distance smaller than 4pc from analytical es- 
timates. The reason for this is the efficient mixing as 
discussed in ij3.2.2l (see also Fig. 

For the survivability of the pre-solar clump we assume 
a quite high density and short lifetime for the Hll-region. 
This leads to a small size, i.e. a small distance for the 
pre-existing clump from the 0-star. The general picture 
would be the ionization front running into the paternal 
filament, from which the 0-star originally formed. On 
a related note, this leads to the assumed high density 
of the ambient surrounding of n = 100 cm~^ (accord- 
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Fig. 8. — The tracer field, which follows the supernova gas for 
case II at t = 14.66 kyr. Color coded is the ratio ^''Al/^'^Al, the 
required value as inferred from CAI measurements is unity (green). 
The length scale is given in units of the radius of the initial cold core 
(i?0 = 0.21 pc). The color coding here is for a 40 Mq progenitor at 
10 pc distance and the supernova gas getting fully mixed within the 
entire supernova bubble. The contours show the density distribu- 
tion. The only regions where the cold gas gets slightly enriched are 
the colliding flanks in the back. Furthermore, no region undergoes 
gravitational collapse. 

ing to the analytical evolution of an Hll-region, see e.g. 
Shu ' 199i[). T h is val ue is well within the observed range. 
iLefloch et al.l ()2002f) estimate t he density in the Trifid 
Nebula to be 50 cm~^ whereas iRubin et all ()2011[ ) find 
values of 87 - 714 cm~^ and 74 - 1041 cm~^ for different 
lines in optical and infrared observations, respectively. In 
turn, the chosen value leads to a high density inside the 
supernova bubble, i.e. a very strong dilution of the en- 
riched material. Thus, optimistic conditions for the sur- 
vival lead to very conservative or pessimistic assumptions 
for the dilution. In addition, 0/B-stars tend to form in 
bigger associations, therefore there would be more ion- 
izing flux than from a single 0-star. Furthermore, we 
neglect the effect of stellar winds. Altogether, the con- 
ditions for the survivability have been chosen quite opti- 
mistically to have a cold core close enough to the 0-star 
to be sufficiently enriched (D « 5pc). Still, as the ob- 
servations around various Hll-regions show, pillar like- 
structures and therefore cold clumps are observed fre- 
quently at these distances. 

On a more technical side, these simulations are still 
only in 2D, so some physical processes might not be 
well resolved. This might apply to the gravitational col- 
lapse of two fronts colliding. However, colliding fronts of 
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these speeds and densities have been shown to become 
gravitationally unstable w ith full 3D simulations (e.g. 
iGritschneder et aLll2009al) . Another short-coming of 2D 
simulations might be the different be havior of instabil- 
ities and mixing in 2D and 3D (e.g. IStone fc Gardinerl 
I2007t ). This should be less of an issue here, as the mix- 
ing appears in the cold phase and is not solely due to 
hydrodynamical instabilities, which would be less well 
represented in 2D (see Fig. The mixing is mainly in- 
fluenced by the cooling, which can be well approximated 
in 2D. Nevertheless, further simulations in full 3D would 
be desirable to address this issues. 

5. CONCLUSIONS & OUTLOOK 

In this study, we focus on the hydrodynamical interac- 
tion of a supernova with the pre-solar cloud core. A re- 
alistic supernova shock wave is assumed to interact with 
the cloud. This leads to the conclusion that it is possible 
to enrich a pre-existing clump sufficiently with ^^Al to 
explain the observed ratio in the Solar System. In or- 
der to survive the previous Hll-region and still be close 
enough to be sufficiently enriched and triggered into col- 
lapse, it turns out that _D « 5pc is an ideal distance. 

The previously stable core gets enriched and is trig- 
gered into gravitational collapse within t « 8kyr after 
the first interaction with the supernova shock-wave. The 
further collapse until the formation of the CAIs, which 
takes places at temperatures below 1800 K, can be esti- 
mated to be shorter than ~ lOkyr. Whether the collapse 
finally goes all the way down to the condensation of the 
CAIs within IBkyr or if there is an intermediate stage 
with a hot disk where the pre-meteorite gas is well mixed 
remains to be investigated. In either case, the CAIs will 
be able to grow fast enough to the observed cm-sizes 
within the background density of a disk, as shown in 
Appendix A. The hot disk, however, can as well only be 
achieved by a very violent (i.e. supernova triggered) col- 
lapse. In addition, the subsequent expansion and cooling 



of the disk will be rapid at these high densities so it is 
preferable if the material is well mixed before. 

On a side note, it is very interesting that the mixing oc- 
curs mainly by material from the supernova shock, which 
already cooled. This cooled portion interacts with the 
pre-existing cold gas in the surrounding of the progen- 
itor star, whereas most of the supernova enriched ma- 
terial stays in the hot component. Therefore, it might 
be possible to determine the efficiency of enrichment by 
supernova feedback in the much larger context of galax- 
ies by investigating the cooling rates within supernova 
shocks. 

Altogether, the precise time and location of the grain 
formation has to be determined by future investigations. 
In this work, we prove that it is possible to trigger a 
pre-existing stable core into collapse by a very nearby 
supernova {D « 5pc). In addition, we show for the first 
time that the collapsing core can be enriched sufficiently 
with SLRs to explain the observed abundances, even if 
the enriched material is not assumed to be in the shock 
initially, but instead distributed homogeneously within 
the entire supernova bubble. 
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APPENDIX 

ASSEMBLY AND GROWTH OF CAIS 

Assuming that Aluminum sticks perfectly to the grains, the growth rate M for a particle of size a at a certain radius 
within the disk with a relative velocity v^ci is 

M = Tra^p^Wrci, (Al) 

where pm is the background density of the CAI-forming material, e.g. ^^Al. If we assume that, as the particle grows 
the density of the particle pp as well as pm remain constant, we can write 

Ana^dpp = na^p^Vrci (^2) 

and solving for d 

a = T^^'ro^ (A3) 

If we ignore the migration of particle, that is, the particle stays in a certain radius r as it grows, the radius a simply 
grows linearly with time 

a = ao + -T^Wroii = ao + kt, (A4) 

where the growth rate k — ^^frci is a function of radius. Since we investigate a longer period of time (t 20kyr), 
the initial particle size is irrelevant, so we can assume ao — 0. Therefore, 

4^"-' = T 



To assemble cm-size objects within 20 kyr, this gives: 



- 1.6 X 10-1^ kms-\ (A6) 



4pp 

Assuming a density pp = lgcm~^ and a relative velocity of Wrci = lkms~^, the resulting Aluminum background 
density has to be p^ = 6.4 x 10~ gcm~^. With the total mass 27 Al of A/27A1 = 7.1 X 10-5 Mq (e.g. iLoddersI 120031) 
this leads to a required total background density of 

Pdisk = 9.0 X 10-13 (A7) 
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which is higher than the density in cloud cores, but can be easily reached in proto-stellar and proto-planetary disks. 

On the one hand this is an optimistic scenario for the growth, since we assume a perfect sticking efficiency. On th e 
other hand, the injected supernova material most likely already contains ^ 0.1/im-sized grains fe.g. lNath et"all[2008h . 
Furthermore, we only take ^''Al into account here, whereas the CAIs are not formed by pure 



